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Perspectives for gapped bilayer graphene polaritonics 
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Bilayer graphene is normally a semimetal with parabolic dispersion, but a tunable bandgap up 
to few hundreds meV can be opened by breaking the symmetry between the layers through an 
external potential. Ab-initio calculations show that the optical response around the bandgap is 
strongly dominated by bound excitons, whose characteristics and selection rules differ from the 
usual excitons found in semiconductor quantum wells. In this work we study the physics of those 
excitons resonantly coupled to a photonic microcavity, assessing the possibility to reach the strong 
and the ultrastrong coupling regimes of light-matter interaction. We discover that both regimes 
are experimentally accessible, thus opening the way for a most promising technological platform, 
combining mid-infrared quantum polaritonics with the tunability and electronic features of graphene 
bilayers. 


I. INTRODUCTION 


Microcavity polaritons are half-light and half-matter 
quasiparticles created by the strong coupling between a 
confined photonic mode and a matter excitation fli- 
Thanks to their hybrid nature they have many unique 
properties, making them the object of intense research 
for applications as different as low threshold lasersJS J 
terahertz emitters quantum simulators 

Since their initial discovery in GaAs quantum wells 
various kinds of polaritonic resonances have been ob¬ 
served in many solid-state cavity quantum electrodynam¬ 
ics setups, ranging from organic microcavities fi3. to in¬ 
tersubband transitions in doped quantum wells [l4| . or 
more recently Landau level systems [l^ . and transition- 
metal dichalcogenides (l^ . 


While microcavity polaritons are not observable in in¬ 
trinsic monolayer graphene, due to the absence of discrete 
excitonic resonances, there have been multiple proposals 
to obtain polaritonic effects in graphene-based systems, 
either by gapping a single graphene layer [ 13 , apply¬ 
ing strong magnetic fields [l8|, or coupling two differ¬ 
ent graphene sheets separated by a dielectric slab di¬ 
ll^. While none of these proposals has been imple¬ 
mented so far, there is little doubt that the observa¬ 
tion of graphene-based microcavity polaritons would be 
an important milestone for present-day nanophotonics. 
Merging two very powerful and versatile technological 
platforms as graphene bilayers and microcavity polari¬ 
tons could lead to major breakthroughs for mid-infrared 
optoelectronics and quantum technologies. In particular 
strong coupling, by increasing the amount of light that 
graphene can absorb or emit, could overcome the main 
obstacle that limits the exploitation of graphene in opto¬ 
electronic devices [25l - l^ . 


In this paper we will investigate the possibility to ob¬ 
serve polaritonic effects in a different graphene-based sys¬ 
tem, that is a bilayer composed of two graphene layers 
in the AB (Bernal) stacking (see the inset in Fig. [I{b)), 
the kind of bilayer obtained by standard graphite exfolia¬ 
tion. As intrinsic graphene, bilayer graphene is normally 
a semimetal, but the interaction between the two layers 


results in parabolic bands around the K and K' points. 
When the symmetry between the two layers is broken, 
e.g., through an externally applied potential, a tunable 
bandgap, up to around 300 meV, can be opened [28l - [^ . 
In Fig. [TJa) we plot the dispersion of the first valence and 
conduction bands for different values of the asymmetry 
potential V. We see that, when a gap is present, the 
shape of the dispersion is not parabolic anymore, but it 
assumes a typical Mexican-hat shape, with the band edge 
on a circle of radius ko in momentum space. Ab-initio 
calculations show that the optical response of gapped bi¬ 
layer graphene (GBG) around the bandgap is strongly 
dominated by bound excitons [i^, Hlj. Those excitons 
are of the Wannier kind, delocalised over many carbon 
atoms on both layers, and they present various peculiar¬ 
ities that distinguish them from excitons in other two 
dimensional materials. In particular, due to its peculiar 
dispersion, GBG presents a van Hove singularity, typical 
of one dimensional systems, in the joint electronic den¬ 
sity of states around the badgap. Not only this diverging 
density of states leads to rather large excitonic features, 
but the Tiw = 2 pseudospin winding number of GBG [l^ 
also strongly modihes optical selection rules [3l|, HI] . 

In the following we will study the physics of those ex¬ 
citons resonantly coupled with the photonic mode of a 
microcavity. Our principal aim will be to estimate the 
maximal achievable ratio between the strength of the 
light-matter interaction, quantified by the vacuum Rabi 
frequency, UIr, and the bare frequency of the excitonic 
resonance, tOx- Such a quantity, that we will call nor¬ 
malised coupling, will give us two fundamental pieces of 
information, allowing us to assess the scientific and tech¬ 
nological interest of this novel solid-state cavity quantum 
electrodynamics setup. The first information is the min¬ 
imal quality factor Q of the photonic resonator which 
allows the system to be in the strong light-matter cou¬ 
pling regime, and thus to sustain polaritonic resonances. 
In order to be in the strong coupling regime the vacuum 
Rabi frequency needs to be larger than the losses, so that 
it becomes possible to spectroscopically resolve the reso¬ 
nant splitting between the two coupled modes [l| . Given 
the large excitonic binding energies of excitons in GBG 
[40l| it is safe to assume that the leading losses are due 
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to the photonic part of the polaritonic field and, remem¬ 
bering that the quality factor of a resonator is the ratio 
between its frequency and its loss rate, the condition to p 
be in the strong coupling regime thus reads g 


> 


4Q' 


CD 



The second information is the extent to which higher or¬ 
der photonic processes are observable in GBG polaritons. 
Applying standard perturbation theory to calculate the 
effects of light-matter coupling we find in fact that the 
Qj-jjgj- perturbative term is proportional to the in¬ 
teraction to the power of n, divided the energy gap to 

( \ n 

— ) . 

^CC J 

The normalised coupling ^ is thus the relevant dimen¬ 
sionless parameter quantifying the coupling between light 
and matter. If such a dimensionless parameter becomes 
non-negligeable, a reg ime referred to as ultrastrong cou¬ 
pling regime [4,11 l44| . many novel physical phenomena 
due to higher order processes become observable, rang¬ 
ing from quantum vacuum radiation and quantum 
phase transitions [H, , to the modification of energy 

transport [3 113 and light emission [sol - fs^ properties, 
to the appearance of cavity assisted chemical and ther¬ 
modynamical effects |53l - l^ . Notice that the transition 
between strong and ultrastrong coupling regimes is a 
smooth crossover, determined more by experimental sen¬ 
sitivity to those novel phenomena that by any intrin¬ 
sic qualitative change in the underlying physics. The 
boundary between the two, that is conventionally fixed 
at ^ =0.1, was in fact the coupling value that allowed 
for tire first observation of ultrastrong coupling effects 
[53 ■ The ultrastrong coupling regime has been since 
then observed in a variety of different systems (b^ - l^ . 
with an actual record of — = 0.87 [OT. Excitons in 
GBG could be very interesting in this regards, a priori 
allowing to reach or even improve such a record. On one 
hand thanks to their peculiar optical properties, leading 
to very strong absorption peaks [40j |. on the other hand 
thanks to the nanometric dimension of the carbon bi¬ 
layer, that can thus be placed very close to the metallic 
surface of metamaterial resonators, achieving very large 
field enhancements [H ill. The rest of this paper is 
structured as follows. In Sec. m we will lay down the 
theory of the coupling between excitons in GBG and the 
photonic mode of a resonator. In Sec. IIIII we will use 
the numerical results of Ref. [13 to estimate the param¬ 
eters we need in order to fully characterise the system. 
With a solid theory and reliable values for the param¬ 
eters, in Sec. m we will then calculate the normalised 
light-matter coupling, showing how both the strong and 
the ultrastrong coupling regimes are experimentally ac¬ 
cessible using GBG. Finally, in Sec. El we will draw 
some conclusions, also pointing out directions for future 
investigations. 



k (in units of V (meV) 


FIG. 1: Left panel: dispersion of the two lowest energy bands 
in graphene bilayer without asymmetry (dash-dotted line), 
an asymmetry parameter V = 200 meV (dashed line) and 
V = 400 meV (solid line). For this last line the Mexican-hat 
shape of the dispersion is evident and the ko point is explicitly 
shown. Right panel: coefficients of the wavefunction as a 
function of the asymmetry, as from Eq. In the inset we 
plot instead a scheme of the bilayer with its tight-binding 
parameters. 


II. THEORY OF LIGHT-MATTER COUPLING 
FOR GBG EXCITONS 


Graphene is an hexagonal lattice of carbon atoms, hav¬ 
ing two atoms per unit cell, usually labelled A and B. 
Around the two inequivalent K and K' points at the cor¬ 
ners of the Brillouin zone it presents a linear dispersion, 
from which most of its exotic electronic properties stem 
[ 13 . A graphene bilayer in AB stacking is obtained by 
putting two graphene monolayers one over the other, such 
that the A atom of one layer is on the top of the B atom of 
the other. In a tight-binding approach the wavefunction 
of a graphene bilayer can be written as a four-component 
vector, describing the wavefunction amplitude in each 
of the two inequivalent sublattices {A, B} in each layer 
{1, 2} [23[33- Around the K point the Hamiltonian can 
be written in the (Ai, Hi, A 2 , H 2 ) basis, keeping only the 
dominant next-neighbour terms, as the direct sum of the 
Hamiltonian for each graphene monolayer (with hopping 
parameter 7 ) plus a term proportional to A, describing 
the interlayer coupling between the dimer sites Hi and 
A 2 
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where U is a parameter describing the asymmetry be¬ 
tween the two layers, cr = UxS-x + o’yCy, with the 
Pauli matrices, cr^ = ± icTy), and Cx^y are linear 

polarisation versors. Such an Hamiltonian can be ana¬ 
lytically diagonalised, leading to the following dispersion 
for the lowest valence and conduction bands as a function 
of the in-plane wavevector k 


= 


+ 2 A 2 


+ - V + A2)72fc2 + ^,(3) 
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whose shape is plotted in Fig. [TJa) for various values of 
V, and whose eigenvectors can be written as 


Fc.k(p) 

F,;,k(p) 


^[al(fc),a2(fc)e*«^a3(A)e*^^a4(fc)e2*«'=],(4) 
^^[-a4(fc), a3(fc)e*®^-a2(fc)e*®^ ai(fc)e^*®''], 


with p the in-plane position vector and S the sample 
surface. In Fig. [Ha) we clearly see the Mexican-hat 
dispersion, with the band edge at a finite value of the 
in-plane wavevector, that can be calculated by differen¬ 
tiating Eq. as 


;,2 _ V*+2V'^A^ 
'^0 — 4 (V‘‘‘+A'-‘) ■ 


( 5 ) 


We expect that the low energy excitons will mainly con¬ 
sist of linear superpositions of electronic transitions lo¬ 
calised around the band edge [t^ , an hypothesys at least 
qualitatively confirmed by the numerical data of Ref. 
[4ll[ . In Fig. (Hb) we thus plot the aj{k) coefficients 
at k = ko, j G {1,2,3,4}. We see that the ai(fco) coef¬ 
ficient is the dominant one, at least for moderate values 
of the asymmetry. In the following, in order to simplify 
the theory, we will thus always assume 

ai(fc) > a2{k), a3{k), a4{k). (6) 

This approximation implies that valence and conduction 
electrons are mainly localised on different layers, in the 
non-dimer sites Ai and B 2 , a conclusion supported by 
the numerical results of Refs. (40l.l4ll|. 

The interaction Hamiltonian can be derived by apply¬ 
ing the usual minimal-coupling substitution p —>■ p + 
eA(r) in Eq. ([2]), where A(r) is the electromagnetic vec¬ 
tor potential and r = {p,z} is the three dimensional 
position vector. We thus obtain 
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the circular polarisation versors and 


( 9 ) 
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where zq is the position of the GBG plane. As firstly 
noted in Ref. [^, the phases of the electron and hole 
wavefunctions in Eg. leading to a pseudospin winding 
number = 2 |42l |. modify selection rules for excitons. 
The sum over the 6k in '!/i(k)e^*®'= implies that only 
excitons with angular dependency proportional to , 
that is having envelope angular momentum ruenv = il, 
are coupled to the crystalline ground state. Under the hy¬ 
pothesis expressed in Eq. m, we can thus consider that 
effectively only the menv = 1 exciton couples around the 
K point. Considering the four-fold spin and valley mul¬ 
tiplicity, and taking care that for excitons around the K' 
point we need to consider the inversion y —>■ —y, leading 
to have the ruenv = ~1 exciton coupled instead [d^, we 
thus obtain the following interaction Hamiltonian 

fl,.. = Mil+ !>,). ( 12 ) 


with 


iFint 



cr 0 \ 

0 cr j ■ 


( 7 ) 


V = *®'=ai(fc)a 3 (fc), (13) 


Note that in Eq. © there is no dipole along z, that 
is consisten with our approximation of considering the 
electrons localised around the non-dimer states, thus not 
allowing for currents normal to the bilayer plane. We 
will thus limit ourselves to the case of normal incidence, 
for which the coupling is maximised as the field lies com¬ 
pletely in the dipole plane. Introducing the second quan¬ 
tised annihilation operators for electrons, Ck, and holes, 
hk, with in-plane wavevector k, the exciton state with 
zero center of mass momentum around the K point can 
be written as 

IV') = (8) 

k 


a parameter depending on the band structure and on the 
exciton wavefunction, and ba the annihilation operator 
for an exciton of polarisation a. 

For the photonic cavity we will use a general model, 
taking into account the possibility to use subwavelength 
confinement in order to increase the coupling. We can 
write the electromagnetic cavity field for a single mode 
of frequency uJc as 

A = J \ (« + at)e, (14) 

where a is the annihilation operator for a photon, e its 
polarisation at the location of the GBG, and the effective 
length, Leff, is defined as [i^ 


with |0) the crystalline ground state. We can now calcu¬ 
late the interaction matrix element in the dipole approx- 


Feff ^/.(zq) 


dz €r{z) u{z)'^, 


(15) 
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FIG. 2: Left panel: fit of the exciton energy as a function 
of the applied potential I4xt. The crosses are numerical data 
taken from the Fig. 4(a) of Ref. and the continuous line 
the numerical fit done using Eq. Cl. Right panel: fit of the 
maximal absorption of the excitonic peak as a function of the 
applied potential I4xt. The crosses are numerical data taken 
from the Fig. 3 of Ref. [^, the continuous line the numerical 
fit done using Eq. CD and Eq. Cl. and the circle is the point 
derived from the integrated absorption in Eq. (I23II . 


where u( 2 ;) is the profile of the electromagentic mode. 
The effective length quantifies to which extent the energy 
of a cavity photon is concentrated around the graphene 
bilayer, and it is thus a crucial element to determine the 
strength of the light-matter interaction: the shorter the 
effective length, the higher the field at the location of the 
excitons, the larger the coupling. Considering only the 
excitonic mode that couples to the photonic one we thus 
finally get 

Hint = , ^ =(a + -I- &), (16) 


from which we can directly read the normalised coupling 
at resonance 


Wa; eo{huixY LeS 


(17) 


III. PARAMETERS ESTIMATION 

In order to extract useful numerical values from 
Eq. (fTTll . able to tell us if either the strong or the ul¬ 
trastrong coupling regimes are observable with GBG, we 
need to have reliable estimates for the coupling parameter 
V, and for the excitonic energy huj^- While it is a priori 
possible to solve the Bethe-Salpeter equation following, 
for example, an approximation scheme similar to the one 
used in Ref. [23, the results obtained in this way are not 
in quantitative agreement with the ab-initio predictions 
from Refs. [40, |4l|. This does not come as a surprise 
given that the method developed in Ref. [t^ assumes 
vanishingly small binding energies, while the binding en¬ 
ergies obtained from ab-initio calculations are sizeable 
fractions of the bandgap. 

In order to get reliable estimates of the achievable nor¬ 
malised couplings, we will thus use data from Ref. [13 to 


obtain fits for hujx and v. Notice that in such a reference 
the results are all presented not as functions of the asym¬ 
metry parameter V that appears in the Hamiltonian in 
Eq. (I2|), but of the physically applied external potential 
14xt- While one can a priori be calculated in function 
of the other, as shown in Ref. [13, resulting in a E of 
the order of 14 xt/ 2 , we do not need to perform such a 
conversion here, because all the parameters we need are 
already given as functions of Eext- For consistency we 
take also the band structure parameters from the same 
reference, even if a few different values are reported in 
the literature [23, that is we will consider 7 = 5.4 eVA 
and A = 0.37 eV. 

The exciton energy (the optical bandgap) is plot¬ 
ted in Fig. 4(a) of Ref. [13 as a function of 14xt- We 
fitted such a plot using a quadratic function 

= C2Eext + ClEext, (18) 

that, as shown in Fig. UKa), gives a very good fit for 
C 2 = —2.9 X 10“^ meV“^ and ci = 0.56. To estimate v 
instead we can exploit the Fig. 3 of Ref. [13, where the 
excitonic absorption for different values of 14 xt is plotted 
assuming an arbitrary HWHM broadening F = 5 meV. 
From Eq. (fT^ we can calculate the absorbed power for a 
normally incident plane wave of frequency u! and vector 
potential A(r) using the Fermi golden rule 


Tabs 


327re^7^A(r)^j^^S' 




(19) 


with p{uj) the density of states of the exciton state, a 
Lorentzian of width F. Dividing Eq. m by the incident 
flux 

Ein = (20) 

we obtain the absorption spectrum 


A{uj) 


647re^7^i^^ 


2567r^a7^z/^ 

JK) 


p{uj), ( 21 ) 


where a is the fine structure constant. At resonance 
p{lS) = and, knowing the excitonic frequencies from 
Eq. (1181) . we can thus fit v against the maximal absorp¬ 
tion for various values of 14xt- In Fig. |2](b) we show that 
a good fit can be obtained with a linear approximation 


= dEext, (22) 

with d = 2.9 X 10^ meV“^m“^. The consistency of such 
a fit can be tested against the value of the integrated ab¬ 
sorption for the lowest bright exciton, IA= 1.24 meV, for 
Eext = 400 meV, explicitly given in Ref. [13 ■ Integrating 
Eq. (I^Tl) we obtain 


4 X lA 


647r^Q;7^^^ 

IKj ’ 


(23) 


where the factor 4 is due to the four-fold exciton degen¬ 
eracy. The obtained value is shown as a circle in Eig. 
Hb). 
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FIG. 3: Left panel: normalised coupling as a function of the 
applied potential I4xt, for a resonant ^ cavity. Right panel: 
normalised coupling as a function of the effective cavity length 
Leff, for I4xt = 400 meV. 


IV. NUMERICAL RESULTS 

Plugging Eq. (118jl and Eq. (122|1 into Eq. m we can 
finally calculate the normalised coupling as a function 
of the external potential. In Fig. ETa) we plot the nor¬ 
malised coupling as a function I4xt for a ^ cavity, that is 
choosing Leg = ^. We see that normalised couplings of 
the order or ^ ~ 0.03 are achievable, that from Eq. ([1]) 
implies that quality factors as low as Q = 10 would 
allow the system to be in the strong coupling regime. 
This is a very bland requirement given that, in the 
mid-infrared, quality factors Q > 10^ have been achieved 
using dielectric cavities based on narrow bandgap IV-VI 
semiconductors [o^. Moreover, graphene layers have 
already successfull y b een coupled to both metallic 
[ 2 ^ and dielectric ^ planar cavities. While those 
experiments were performed at shorter wavelengths, 
they are nevertheless important to assess the feasibility 
of such a kind of device. In Fig. [IJb) we plot instead the 
normalised coupling as a function of the effective length 
Leff, normalised to the free space wavelength A. We 
see that the ultrastrong coupling regime is achievable 
already for subwavelength confinements Leff — 0.03A, 
and that ver y h igh values are achievable for tighter ones. 
From Ref. [6^, where a simple condenser-type cavity 
is studied, we see that Leff small enough are achievable 
while still having a sizeable quality factor. Notice that 
linear subwavelength confinements, of the order of 10“^ 
in the mid-infrared, corresponding to ^ ~ 0.2, are 
commonly achieved in metallic structures or, more 
rece ntly , in localised surface phonon polariton samples 
[ 711 1721 , where such low mode volumes are accompanied 
by quality factors Q > 100, over an order of magnitude 
larger of what necessary to resolve the polaritonic 
resonances. Mid-infrared subwavelength resonators have 
also already been coupled to single graphene layers (t^ . 
Notice that, contrary to other kinds of mid-infrared 
polaritonic structures based on intersubband transitions, 
where the effective field confinement is limited by the 
fact that it is not possible to put the two dimensional 
electron gas closer than a few tens of nanometers from 
the metallic plasmonic surface [b^, in our case the GBG 


can effectively be placed in contact with it, at the same 
time increasing the coupling and allowing to use the 
metallic cavity as one of the gates needed to create the 
bandgap. 


V. CONCLUSIONS AND PERSPECTIVES 

We studied the physics of excitons in GBG resonantly 
coupled to a photonic inicrocavity, assessing the possi¬ 
bility to reach the strong and the ultrastrong regimes of 
cavity quantum electrodynamics, and we discovered that 
both of them are within experimental reach. We hope 
that such a result will foster further interest in the study 
of light-matter coupling in GBG, eventually leading to 
the realisation of graphene-based polaritonic microcav¬ 
ities. Such devices could empower novel technological 
breakthroughs in the mid-infrared region of the electro¬ 
magnetic spectrum, by joining the unique physical prop¬ 
erties of tunable gapped graphene and strongly coupled 
polaritons. 

While the theory presented here is mainly phenomeno¬ 
logical, further studies, coadiuvated by ab-initio calcula¬ 
tions, could give a clearer picture of the unique properties 
of those polaritons. In particular it is important to notice 
that the absolute value of the coupling Ll/j can become 
comparable with the excitonic binding energy, thus oblig¬ 
ing a consistent theory to consider also the coupling with 
higher lying excitonic modes and with the continuum of 
electron states. This kind of effects, meviously studied in 
quantum well-based microcavities [74l476| . amounts to a 
cavity induced modification of the exciton wavefunction, 
and it leads to an alteration of the spectral features of 
the polaritonic resonances. 

It is worthwhile to notice that the theory we devel¬ 
oped in Sec. lUis based on an effective minimal coupling, 
inspired to what usually done in the theory of micro¬ 
cavity polaritons in quantum wells [ 3 , 123 . Still, a dia¬ 
tribe rece ntly arose on the pertinence of such a model for 
graphene whose crucial point is the presence, or 

lack of thereof, of the diamagnetic A(r)^ term in the min¬ 
imal coupling Hamiltonian. While the presence of such 
a term does not change the intensity of the coupling, 
and it does not thus modify the results of the present 
work, it still has a non-negligible impact on the system. 
Not only it shifts the polaritonic energies [13, but for 
high normalised couplings its presence causes an effec¬ 
tive decoupling between light and matter [ 10 , [13 , while 
its absence leads to an instability of the ground state and 
a subsequent superradiant phase transition [10, EJ [lOl . 
While the solution of this diatribe lies well outside the 
scope of the present paper, we wish to point out that 
the very large couplings achievable in GBG polaritons, 
and the consequent very different predictions for observ¬ 
ables with and without the A(r)^ term, would allow a 
definitive solution of this conundrum. 
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